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Abstract 

We investigate the possibility of detecting solar electron antineutrinos with the 
KamLAND experiment. These electron antineutrinos are predicted by spin-flavor 
oscillations at a significant rate even if this mechanism is not the leading solution to 
the SNP. 

KamLAND is sensitive to antineutrinos originated from solar neutrinos. From 
KamLAND negative results after 145 days of data taking, we obtain model inde- 
pendent limits on the total flux of solar electron antineutrinos < 1.1 — 3.5 x 
10^ cm~^ s~^, more than one order of magnitude smaller than existing limits, and 
on their appearance probability P < 0.15% (95% CL). 

Assuming a concrete model for antineutrino production by spin-flavor precession, 
this upper bound implies an upper limit on the product of the intrinsic neutrino 
magnetic moment and the value of the solar magnetic field fiB < 2.3 x 10~^^ MeV 
95% CL (for LMA {Am^,tan^e) values). 

Limits on neutrino transition moments are also obtained. For realistic values 
of other astrophysical solar parameters these upper limits would imply that the 
neutrino magnetic moment is constrained to be, in the most conservative case, 
/x<3.9 X 10-12 (95% QL) foj. ^ relatively smah field B = 50 kG. For higher 
values of the magnetic field we obtain: /x<9.0 x 10"^^ fiB for field B = 200 kG and 
//<2.0 X 10-13 /is for field B = 1000 kG at the same statistical significance. 
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1 Introduction 



Evidence of eelctron antineutrino disappearance in a beam of antineutrinos in the 
KamLAND experiment has been recently presented [1] . The analysis of these results 
[1,2] in terms of neutrino oscillations have largely improved our knowledge of neu- 
trino mixing in the LMA region. The results appear to confirm in a independent way 
that the observed deficit of solar neutrinos is indeed due to neutrino oscillations. The 
ability to measure the LMA solution, the one preferred by the solar neutrino data at 
present, "in the lab" puts KamLAND in a pioneering situation: after these results 
there should remain little doubt of the physical reality of neutrino mass and oscilla- 
tions. Once neutrino mass is observed, neutrino magnetic moments are an inevitable 
consequence in the Standard Model and beyond. Magnetic moment interactions arise 
in any renormalizable gauge theory only as finite radiative corrections: the diagrams 
which contribute to the neutrino mass will inevitably generate a magnetic moment 
once the external photon line is added. 

The KamLAND experiment is the successor of previous reactor experiments ( 
CHOOZ [3], Palo Verde [4]) at a much larger scale in terms of baseline distance and 
total incident fiux. This experiment relies upon a 1 kton liquid scintillator detector 
located at the old, enlarged, Kamiokande site. It searches for the oscillation of 
antineutrinos emitted by several nuclear power plants in Japan. The nearby 16 
(of a total of 51) nuclear power stations deliver a Pg Aux of 1.3 x lQ^cm~^s~^ for 
neutrino energies E,y > 1.8 MeV at the detector position. About 85% of this flux 
comes from reactors forming a well defined baseline of 139-344 km. Thus, the flight 
range is limited in spite of using several reactors, because of this fact the sensitivity 
of KamLAND increases by nearly two orders of magnitude compared to previous 
reactor experiments. 

Beyond reactor neutrino measurements, the secondary physics program of Kam- 
LAND includes diverse objectives as the measurement of geoneutrino flux emitted by 
the radioactivity of the earth's crust and mantle, the detection of antineutrino bursts 
from galactic supernova and, after extensive improvement of the detection sensitivity, 
the detection of low energy ^Be neutrinos using neutrino-electron elastic scattering. 

Moreover, the KamLAND experiment is capable of detecting potential electron 
antineutrinos produced on fly from solar neutrinos [5]. These antineutrinos are 
predicted by spin-flavor oscillations at a significant rate if the neutrino is a Majorana 
particle and if its magnetic moment is high enough [6, 7]. In Ref.[5] as been remarked 
that the flux of reactor antineutrinos at the Kamiokande site is comparable, and in 
fact smaller, to the flux of ^B neutrinos emitted by the sun ,^{^B) ~ 5.6x10^ cm s 
[1, 8, 9]. Their energy spectrum is important at energies 2—4 MeV while solar neutrino 
spectrum peaks at around 9—10 MeV. As the inverse beta decay reaction cross section 
increases as the square of the energy, we would expect nearly 10 times more solar 
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electron antineutrino events even if the initial fluxes were equal in magnitude. 

The publication of the SNO results [9, 10] has already made an important break- 
through towards the solution of the long standing solar neutrino [11-14] problem 
(SNP) possible. These results provide the strongest evidence so far (at least until 
KamLAND improves its statistics) for flavor oscillation in the neutral lepton sector. 

The existing bounds on solar electron antineutrinos are strict. The present upper 
limit on the absolute flux of solar antineutrinos originated from neutrinos is 
[5, 16, 17] <&i7(^i?) < 1.8 X 10^ cm^'^ s^^ which is equivalent to an averaged conversion 
probability bound of P < 3.5% (SSM-BP98 model). There are also bounds on their 
differential energy spectrum [16]: the conversion probability is smaller than 8% for 
all Ee^vis > 6.5 MeV going down the 5% level above £?e,ms — 10 MeV. 

The main aim of this work is to study the implications of the recent KamLAND 
results on the determination of the solar electron antineutrino appearance probability, 
independently from concrete models on antineutrino production. The structure of 
this work is the following. In section 2 we discuss the main features of KamLAND 
experiment that are relevant for our analysis: The salient aspects of the procedure we 
are adopting and the results of our analysis are presented and discussed in sections 
3. In Section 4 we apply the results we obtained in a particular model for the solar 
magnetic field, we obtain bounds on the values of the intrinsic neutrino transition 
magnetic moments. Finally, in section 5 we draw our conclusions and discuss possible 
future scenarios. 



2 The computation of the expected signals 
2.1 The KamLAND signal 

Electron antineutrinos from any source, nuclear reactors or solar origin, with energies 
above 1.8 MeV are measured in KamLAND by detecting the inverse /3-dccay reaction 
Ve + p ^ n + e'^. The time coincidence, the space correlation and the energy balance 
between the positron signal and the 2.2 MeV 7-ray produced by the capture of a 
already-thermalized neutron on a free proton make it possible to identify this reaction 
unambiguously, even in the presence of a rather large background. 

The two principal ingredients in the calculation of the expected signal in Kam- 
LAND are the corresponding flux and the electron antineutrino cross section on pro- 
tons. The average number of positrons Ni originated from the solar source which are 
detected per visible energy bin AEi is given by the convolution of different quantities: 

r f'CO f'CO 

Ni = Qo dEe dEle{Ee)R{Ee,El) dE,jp{E,j)^E,j)a{E,j, E^) (1) 

J AEi Jo J El 

where Qq is a normalization constant accounting for the fiducial volume and live 
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time of the experiment, p. Expressions for the electron antineutrino capture cross 
section a{ETj,El) are taken from the hterature [18, 19]. The matrix element for this 
cross section can be written in terms of the neutron half-life, we have used the latest 
published value ti/2 = 613.9 ± 0.55 [17]. The functions e{Ee) and R(Ee,El) are 
the detection efficiency and the energy resolution function. We use in our analysis 
the following expression for the energy resolution in the prompt positron detection 
(j{E(,) = 0.0062 + 0.065-y-Ee . This expression is obtained from the raw calibration 
data presented in Ref.[20]. Note that we prefer to use this expression instead of the 
much less accurate one given in Ref.[l]. Moreover, we assume a 408 ton fiducial 
mass and the detection efficiency is taken independent of the energy [1], e = 80%. 
In order to obtain concrete limits, a model should be taken which predict p and its 
dependence with the energy. For our purpose it will suffice to suppose p a constant. 
This is justified at least in two ) if the energy range /S.E over which we perfom 

the integration is small enough so the variation of the probability is not very large, 
or b) if we reinterpret p as an energy-averaged probability, note that, in a general 
case, this is always true because the un-avoidable convolution with a finite energy 
resolution, (see Expression 10 in Ref.[7]): 

p^E = I dE a{E)^{E)PT,{E)/ j dE a{E)^{E). (2) 

J AE JAE 

Let us finally note that independently of the reasons above, upper limits to be ob- 
tained on continuation are still valid even if the antineutrino probabilities are signifi- 
cantly different from constant: if we take p = maxAc PwiE) the expected antineutrino 
signal computed with p will be always larger than the signal obtained inserting the 
full probability. 

The results of our simulation are summarized in Fig. 5 where we show the "solar" 
positron spectrum obtained assuming the shape of the neutrino flux and a total 
normalization 10~^ x $(^5) which means an overall Ug — Vg conversion probability 
P ~ 1%. 

In addition we have computed the expected signal coming from antineutrino reac- 
tors. A number of short baseline experiments (See Rcf.[21, ?] and references therein) 
have previously measured the energy spectrum of reactors at distances where oscilla- 
tory effects have been shown to be inexistent. They have shown that the theoretical 
neutrino flux predictions are reliable within 2% [22] . The effective flux of antineutri- 
nos released by the nuclear plants is a rather well understood function of the thermal 
power of the reactor and the amount of thermal power emitted during the fission 
of a given nucleus, which gives the total amount, and the isotopic composition of 
the reactor fuel which gives the spectral shape. Detailed tables for these magni- 
tudes can be found in Ref. [21, ?]. For a given isotope the energy spectrum can be 
parametrized by an exponential expression [18] where the coefficients depend on the 
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nature of the fissionable isotope (see Ref.[21] for expHcit values). Along the year, 
between periods of refueling, the total effective flux changes with time as the fuel is 
expended and the isotope relative composition varies. We take the average of the 
relative fission yields over the live time as given by the experiment: 235jj _ 57%^ 
238[/ = 7 8%^ 239 p„ ^ 30%^ 241 ^ 5 7% ^^.^^^ obtain the expected number 
of events at KamLAND, we sum the expectations for all the relevant reactor sources 
weighting each source by its power and distance to the detector (table II in Ref. [21]), 
assuming the same spectrum originated from each reactor. Wc sum over the nearby 
power reactors, we neglect farther Japanese and Korean reactors and even farther 
rest-of-the-world reactors which give only a minor additional contribution. 



We will obtain upper bounds on the solar electron antineutrino appearance proba- 
bility analyzing the observed KamLAND rates in three different ways. In the first 
one, we will make a standard analysis of the observed and expected solar signals 
in the 13 prompt positron energy bins considered by KamLAND [1]. In the second 
and third cases we will apply Gaussian and poissonian probabilistic considerations 
to the global rate seen by the experiment and to the individual event content in the 
highest energy bins (Ee > 6 McV) where KamLAND observes zero events. This null 
signal makes particularly simple the extraction of statistical conclusions in this case. 

3.1 Analysis of the KamLAND Energy Spectrum 

Here we fully use the binned KamLAND signal (see Fig. 5 in Ref.[l]) for estimat- 
ing the parameters of solar electron antineutrino production from the method of 
maximum- likelihood. We minimize the quantity 



where the first term correspond to the contribution of the first nine bins where the 
signal is large enough and the use of the Gaussian approximation is justified. The 
second term correspond to the latest bins where the observed and expected signals 
are very small and poissonian statistics is needed. The exphcit expressions are: 



3 Analysis and Results 



X — Xi=l,9 + Xi=10,13 



(3) 




(4) 



j=l,9 




= 2 
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The quantities Si are the observed bin contents from KamLAND. The theoretical 
signals are in principle a function of three different parameters: the solar electron an- 
tineutrino appearance probability p and the neutrino oscillation parameters {Am?, 6). 
Both contributions, the contribution from solar antineutrinos and that one from solar 
reactors, can be treated as different summands: 

Si{p,Am\e) = 5f'""(p) + 5r"^*°''(Am^^). (6) 

According to our model, the solar antineutrino appearance probability p is taken as 
a constant and we can finally write: 

Si{p,Am'^,9) = px 5° + 5^"^*''^(Am^^). (7) 

In this work we will take for the minimization values of the oscillation parameters 
those obtained when ignoring any solar antineutrinos (LMA solution Am^ = 6.9 x 
10~^ eV^, sin^ 6 = 1 from Ref.[l]) and we will perform a one-parameter minimization 
with respect p. This approximation is well justified because the solar antineutrino 
probability is clearly very small. We avoid in this way the simultaneous minimization 
with respect to the three parameters {p, Am? , 9) . 

We perform a minimization of the one dimensional function x^ip)- to test a 
particular oscillation hypothesis against the parameters of the best fit and obtain the 
allowed interval in p parameter space taking into account the asymptotic properties 
of the likelihood function, i.e. log JC — log JCmin behaves asymptotically as a with 
one degree of freedom. In our case, the minimization can be performed analytically 
because of the simple, lineal, dependence. A given point in the confidence interval is 
allowed if the globally subtracted quantity fulfills the condition Ax^ = (p) ~ Xmm < 
xl{CL). Where Xn=i(90%' 9^%, ...) = 2.70,3.84,.. are the quantiles for one degree 
of freedom. 

Restricting to physical values of p, the minimum of the x^ function is obtained 
for p = Q. The corresponding confidence intervals are p < 4.5% (90% CL) and 
p < 7.0% (95% CL). We have explicitly checked, varying the concrete place where 
the division between "Gaussian" and "poissonian" bins is established in Expression 3, 
that the values of these upper limits are largely insensitive to details of our analysis. 
In particular, similar upper limits are obtained in the extreme cases: if Gaussian or 
poissonian statistics is employed for all 13 bins. These upper limits are considerably 
weaker than those obtained in the next section. One possible reason for that is that 
they are obtained applying asymptotic general arguments to the x^ distribution, 
stronger, or more precise limits could be obtained if a Monte Carlo simulation of the 
distribution of the finite sample distribution is performed (where the boundary 
condition p > should be properly included) . 



6 



3.2 Analysis of the global rate and highest energy bins 

We can make an estimation of the upper bound on the appearance solar electron 
antineutrino probability simply counting the number of observed events and sub- 
tracting the number of events expected from the best-fit oscillation solution. For our 
purposes this difference, which in this case is positive, can be interpreted as a hypo- 
thetical signal coming from solar antineutrinos (Amg = 6.9 x 10^ ey^,sin^^o = 1)- 

Ssolar = PX Ssoig_r = Sobs - Sreacti^Ull, sin^ 6*0). (8) 

Putting [1] Sobs = 54.3 ± 7.5 and S'react ( Amg, sin^ 6'o) = 49 ± 1.3, we obtain Sgbs — 
Sreact < 64.8 (67.2) at 90 (95)% CL. Prom these numbers, the corresponding limits 
on solar electron antineutrino appearance probability are p < 0.45%, 0.52% at 90 
or 95% CL. These limits are valid for the neutrino energy range E,y ^ 2 — 8 MeV. 
In this case, due to the large range, the limits are better interpreted as limits on an 
energy-averaged probability according to expression 2. 

In a similar approach, we use on continuation the binned KamLAND signal corre- 
sponding to the four highest energy bins (see Fig. 5) which, as we will see, provide the 
strongest statistical significance and bounds. The reason for that is that the exper- 
iment KamLAND does not observe any signal here and, furthermore, the expected 
signal from oscillating neutrinos with LMA parameters is negligibly small. 

Due to the small sample, we apply Poisson statistics to any of these bins and use 
the fact that a sum of Poisson variables of mean /Xj is itself a Poisson variable of mean 
^ . The background (here the reactor antineutrinos) and the signal (solar electron 
antineutrinos) are assumed to be independent Poisson Random variables with known 
means. 

If no events are observed, and, in particular, no background is observed, the 
unified intervals [17,23] [0,ecL] are [0,2.44] at 90% CL and [0,3.09] at 95% CL. 
Prom here, we obtain p x S'q"'"'' < ecL or p < ech/ S^"^"''^ . Taking the expected 
number of events in the first 145 days of data taking and in this energy range (6-8 
MeV) we obtain: p < 0.12% (90% CL) and p < 0.15% (95% CL). 

4 A model for solar antineutrino production 

The combined action of spin flavor precession in a magnetic field and ordinary neu- 
trino matter oscillations can produce an observable flux of V^rs from the Sun in the 
case of the neutrino being a Majorana particle. In the simplest model, where a thin 
layer of highly chaotic of magnetic field is assumed at the bottom of the convective 
zone (situated at i? ~ O.?/?©), the antineutrino appearance probability at the exit 
of the layer Piv) is basically equal to the appearance probability of antineutrinos at 
the earth [6, 7] ( see also Refs.[24] for some recent studies on RSFP solutions to the 
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Solar Neutrino Problem). The quantity PiT') is in general a function of the neutrino 
oscillation parameters {Am?, 9), the neutrino intrinsic magnetic moment and also of 
the neutrino energy and the characteristics and magnitude of the solar magnetic field. 
However, in a accurate enough approximation, such probability can be factorized in 
a term depending only on the oscillation parameters and another one depending only 
on the spin-flavor precession parameters: 



where Pg^ is the e — ji solar conversion probability. We will assume in this work 
the LMA central values for {Arti^^O) obtained from recent KamLAND data and 
which are compatible with the SNO observations in solar neutrinos [25], we will take 
Pen{Am?,9) ~ (Pen)^^^'^'^'-^ c^; 0.4. The second factor appearing in the expression 
contains the effect of the magnetic field. This quantity depends on the layer width 
Ar (~ O.li?0) and = |Lo^^(i?^), where {B'^) the r.m.s strength of the magnetic 
field and Lq is a scale length {Lq ~ 1000 km). For small values of the argument 
we have the following approximate expression which is accurate enough for many 
applications 



the solar astrophysical factor k = 2/3Peij,LoAr is numerically k^^^ ~ 2.8 x 10~' 



MeV"^. Upper limits on the antineutrino appearance probability can be translated 
into upper limits on the neutrino transition magnetic moment and the magnitude of 
the magnetic field in the solar interior. The results of the Formula 9 can be seen in 
Figure 2. An upper bound p < 0.15 — 0.20% (95% CL) implies an upper limit on the 
product of the intrinsic neutrino magnetic moment and the value of the convective 
solar magnetic field as fiB < 2.3 x 10~^^ MeV (95% CL). In Fig.2 we show the 
antineutrino probability as a function of the magnetic moment fi for fixed values of 
the magnitude of the magnetic field. For realistic values of other astrophysical solar 
parameters {Lq ~ 1000 km, Ar ~ 0.1 Rq), these upper limits would imply that 
the neutrino magnetic moment is constrained to be, in the most desfavourable case, 
H<3.9 X 10~^^ hb (95% CL) for a relatively small field i? = 50 kG. Stronger limits are 
obtained for slightly higher values of the magnetic field: //<9.0 x 10~^^ ij,b (95% CL) 
for field B = 200 kG and /i<2.0 x IQ-^^ (95% ^j. ggj^j 5 ^ ^qOO kG. Let us 
note that these assTimcd vahics for the magnetic field at the base the solar convective 
zone are relatively mild and well within present astrophysical expectatives. 



In summary in this work we investigate the possibility of detecting solar antineutrinos 
with the KamLAND experiment. These antineutrinos are predicted by spin-flavor 




■)] 



(9) 
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Conclusions 
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solutions to the solar neutrino problem. 

The KamLAND experiment is sensitive to potential antineutrinos originated from 
solar *B neutrinos. We find that the results of the KamLAND experiment put rel- 
atively strict limits on the flux of solar electron antineutrinos < 1.1 — 3.5 x 
10^ cm~^ s~^, and their energy averaged appearance probability (P < 0.15 — 0.50%). 
These limits are largely independent from any model on the solar magnetic field 
or any other astrophysical properties. As we remarked in Section 2.1, these upper 
limits on antineutrino probabilities and fluxes arc still valid even if the antineutrino 
probabilities are significantly different from constant. 

Next we assume a concrete model for antineutrino production where they are 
produced by spin-flavor precession in the convective solar magnetic field. In this 
model, the antineutrino appearance probability is given by a simple expression as 
P{17) = K fi^{B'^) with K^^^^ ~ 2.8 x 10""^^ MeV"^ In the context of this model 
and assuming LMA central values for neutrino oscillation parameters (Am^ = 6.9 x 
10~^ eV^, sin^ ^ = 1) [1], the upper bound p < 0.15% (95% CL) implies an upper 
limit on the product of the intrinsic neutrino magnetic moment and the value of the 
convective solar magnetic field as n B < 2.3 x 10~^^ MeV (95% CL). For realistic 
values of other astrophysical solar parameters these upper limits would imply that 
the neutrino magnetic moment is constrained to be, in the most desfavourable case, 
/i<3.9 X 10"^2 (95% for ^ relatively small field B = 50 kG. For slightly 
higher values of the magnetic field: fi<9.0 x 10"^^ i^b (95% CL) for field B = 200 
kG and /x<2.0 x 10"^^ hb (95% CL) for field B = 1000 kG. These assumed values 
for the magnetic field at the base the solar convective zone are relatively mild and 
well within present astrophysical expectatives. 
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Figure 1: The KamLAND positron spectra from reactor antineutrinos (from Fig. 5 in 
Ref. [1]): measured (145.5 days), MC expectations in absence of oscillations and best fit 
including neutrino oscillations (Am^ = 6.9 x 10~^ eV^, sin^ ^ = 1, respectively points with 
error-bars, triangles and stars). The "solar" positron spectrum (black solid squares) ob- 
tained assuming the shape of the neutrino flux and a total normalization 10^^ x $(^i?) 
(that is, an overall Ue — z/g conversion probability P ~ 1%). 
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Figure 2: The solar antineutrino appearance probability p as a function of the transition 
neutrino magnetic moment, in units of Bohr magnetons /z^, for fixed values of the r.m.s 
solar magnetic field (Formula 9). From left (solid) to right (dashed), curves correspond 
to B = 1000,200,50 kG. From the curves, an upper limit p < 0.15% implies /i < 1.9 x 
10~^^/iB,9.0 X 10~^^/U5,3.0 X lO^^^yU^ respectively for each of the magnetic field above. 
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